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ABSTRACT 


There exists a widely spread notion that gravitational effects can strongly violate 
global symmetries. If this is correct, it may lead to many important consequences. 
We will argue, in particular, that nonperturbative gravitational effects in the axion 
theory lead to a strong violation of CP invariance unless they are suppressed by an 
extremely small factor g < 10“®^. 

One could hope that this problem disappears if one represents the global symme¬ 
try of a pseudoscalar axion held as a gauge symmetry of the Ogievetsky-Polubarinov- 
Kalb-Ramond antisymmetric tensor held. We will show, however, that this gauge 
symmetry does not protect the axion mass from quantum corrections. The amplitude 
of gravitational ehects violating global symmetries could be strongly suppressed by 
e~^, where S is the action of a wormhole which may eat the global charge. Unfortu¬ 
nately, in a wide variety of theories based on the Einstein theory of gravity the action 
appears to be fairly small, S = 0(10). 

However, we have found that the existence of wormholes and the value of their 
action are extremely sensitive to the structure of space on the nearly Planckian scale. 
We consider several examples (Kaluza-Klein theory, conformal anomaly, terms) 
which show that modihcations of the Einstein theory on the length scale I < lOMp^ 
may strongly suppress violation of global symmetries. We have found also that in 
string theory there exists an additional suppression of topology change by the factor 

e . This effect is strong enough to save the axion theory for the natural values of 
the stringy gauge coupling constant. 
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1 Introduction 


The most elegant way to solve the strong CP violation problem is given by the Peccei-Qninn 
mechanism |]l|. This mechanism is based on the assnmption that there exists a complex scalar 
held <h(a:) = ^ which after spontaneons symmetry breaking can be represented as 

exp(^^^^^. The Goldstone held a{x) (axion) has the conpling similar to the 

famons 6'-term . Nonpertnrbative ehects in QCD lead to appearance of the condensate 

{F^yF^y), and to the ehective potential of the axion held proportional to Aq^^( 1 — cos(0 + ^)). 


This potential has a minimnm at ^ = —9. In this minimnm the terms 


327r^/o 


FyyF^^'^ and 


y^F^yF^'" cancel each other, and strong CP violation disappears. This ehect gives the axion a 


small mass rUa 


^QCD 

fo ■ 


In addition to providing a possible solntion to the strong CP violation problem, invisible axion 
held PI is one of the best dark matter candidates p. It natnrally appears in all phenomenological 
models based on snperstring theory 0]. Axion held possesses many interesting properties near 
black holes p. Finally, axions may be responsible for a possible existence of wormholes in 
the baby-nniverse theory |^. Therefore there exists an extensive literatnre on axions. This 
literatnre inclndes at least two diherent formnlations of the axion theory, which are not completely 
eqnivalent, and several modihcations of these formnlations (for a review see P). 

The axion theory has many problems. First of all, it is not easy to make this theory compatible 
with cosmology. If the spontaneons symmetry breaking towards a state with fo ^ 0 occnrs 
after the end of inhation, then the standard axion model is compatible with cosmological and 
astrophysical constraints only if 10^° GeV < /o ^ 10^^ GeV |^. Recent investigation with an 
acconnt of cosmological ehects of the axion strings snggests that the npper bonnd may be even 
more tight, so that the “axion window” becomes almost closed, 10^*^ GeV < /o < IQii GeV 
0. On the other hand, if the spontaneous symmetry breaking occurs during inhation, then the 
constraint fo < 10^^ disappears |]Tp, but typically it implies that the Hubble constant at the end 
of inhation should be sufficiently small, FT < 10® GeV |^. Inhationary models of this type 
can be easily suggested |Tp, but one should keep in mind that not every inhationary model 
satishes this condition. 


Three years ago it was pointed out [jlp that the axion theory faces another difficult problem, 
which we are going to discuss in this paper. 
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The standard potential in the axion theory (ignoring small QCD corrections) is given by 

l/o($) = A(|$p-/oV2)^ (1) 


In this approximation the axion is massless dne to the global symmetry $ —»• However, 

there are some reasons to expect that nonpertnrbative qnantnm gravity effects do not respect 
global symmetries. The simplest way to nnderstand it is to remember that global charges can be 
eaten by black holes, which snbseqnently may evaporate. One may expect that a similar effect 
can occnr becanse of a nonpertnrbative formation and evaporation of “virtnal black holes” in 
the presence of a global charge. A somewhat more developed (even thongh still very specnlative) 
approach is based on investigation of wormholes, which may take a global charge from onr nniverse 
to some other one. Indeed, it was claimed in p, |^, [Tp] that snch effects do actnally take place, 
and can be described by additional terms (vertex operators) in the effective Lagrangian which 
break the global symmetry. 


As an example, one may consider the terms of the following type 

12m, ^4—2m,H-n 


V);($) = 9n 


$ 




+ h.c. 


( 2 ) 


where gn is some dimensionless constant. Naively, one conld expect that these operators shonld 
be at least of the hfth order in <h, so that they shonld be snppressed by Mp in the denominator, 
with n > 0. The anthors of concentrated on the simplest (and the most dangerons) term 
They have shown that, for /o ~ 10^^ GeV, this term destroys the standard solntion 


95- 


of the strong CP problem. Indeed, this term changes the shape of the effective potential and 


-6. If ^ changes by more than 10 the corresponding 


moves its minimnm away from ^ = 
effects of CP violation become too strong. In order to avoid snch effects, one shonld have an 
extremely small conpling constant g of the symmetry breaking operator g^ ^^^ g^ < 

for /„ ~ 1012 GeV H. Thns, instead of the problem of explaining why the angle 6 in the theory 
of strong interactions is smaller than 10“®, we mnst explain now why some other parameter is 
smaller than This does not look like a fair trade! 


In fact, the sitnation is even more complicated. The idea to consider only the terms containing 
Mp in the denominator was based on the assnmption that the qnantnm gravity effects shonld 
be snppressed in the limit Mp —>■ cxo. Indeed, the n-loop qnantnm gravity corrections contain 
factors Mp"'. However, the effects we are interested in are nonperturbative. These effects may 
give rise to vertex operators of the type of giMp{^ + ^*), or other operators which do not contain 
Mp in the denominator |^, The way to see it is, e.g., to consider the averages of the type 
Mp((<f) + $*)). One can show that in the presence of wormholes snch terms do not vanish, bnt 
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they are suppressed by the same exponential factor e as the terms [0- Here S 


is the action of a wormhole which can “eat” a unit of a global charge associated with the field 


<h. This implies that if the effective vertex operators 


Mp 


appear in the theory due to 


nonperturbative gravitational effects, one may expect that the operators + $*) should 

appear as well, with a comparable coupling constant, ~ (75 ~ e~^. 


The vertex operator giM^{^ + $*) is most dangerous for the axion physics. One can easily 
show, by analogy with |]^, that this term leads to a strong CP violation unless gi < 10~^^ iQi^ceV 
This constraint is almost thirty orders of magnitude stronger that the constraint following from 
the investigation of the operators ^ ■ For comparison of our constraint with the results 

of our future calculations of the wormhole action it is convenient to express this constraint (for 
/o ~ 10^^ GeV) in the form gi < 


Note, that this constraint depends on the value of /o, but not too strongly: it is proportional 
to /o- Thus, for /o ~ 10^° GeV one should have gi < whereas for /q ~ 10^® GeV our 

constraint is gi < In what follows we will suppose, for definiteness, that /o = 10^^ GeV, 

even though, as we have already emphasized, /o in the axion theory may be either two orders of 
magnitude smaller, or much greater than 10^^ GeV. 


Similar arguments are valid for other theories possessing global symmetries. For example, 
it was shown in that the theory of cosmic textures may work only if the constant g in 
the term ( 75 is extremely small; g^ < 10“®^. One can easily show that this constraint 
becomes even much stronger if one takes into account the above-mentioned terms linear in $: 
gi < 10-103 ^ e-23T 


The same situation appears in the so-called “natural inflation” model [|T^. In this model it 
is assumed that the effective potential has the form (j^) with /o > Mp, and then the Goldstone 
field, just like the axion field, acquires mass h ?//o ~ GeV. This can be achieved in a natural 

way, e.g., for A ~ lO^o GeV and /o ~ Mp. However, the term 5 fiMp(<F-|- <F*) will destroy this nice 
picture unless the coupling gi is extremely small, gi < IQ-®. As compared with the constraints on 
the models of axions and textures, this condition looks relatively mild. Still, the existence of this 
additional constraint is rather disappointing. (On the other hand, it would be quite encouraging 
to hnd a natural mechanism which would lead to the gravity-induced terms 5 fiM|($-F$*) with 
9i ^ 10 “®, since some terms breaking the global symmetry are necessary in this scenario.) 


As we already mentioned, the main reason why quantum gravity may break global symmetries 
is associated with the possibility that the global charge can be eaten by wormholes (or virtual 
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black holes) and taken away from onr nniverse. It is commonly believed, however, that local 
charges, snch as electric or magnetic one, cannot disappear, and therefore qnantnm gravity does 
not break local symmetries. The reason can be easily nnderstood if one think abont electric (or 
magnetic) charges falling into black holes. Becanse of the Ganss law, the flnx of electric held 
cannot disappear when the charge falls into a black hole. Charged black holes cannot evaporate 
entirely and take the electric charge away from onr nniverse. Instead of that, they eventnally 
form charged extreme black holes which do not evaporate any fnrther. 


It wonid be very tempting to nse a similar mechanism to save the axion theory. Indeed, it is 
well known that the theory of a massless psendoscalar axion held in a topologically trivial space 
is eqnivalent to the theory of an antisymmetric tensor gange held 6 ^ 1 ,. This held, which was 
introdnced by Ogievetsky and Polnbarinov and later also by Kalb and Ramond [^, natnrally 
appears in string theory. It is related to the held a by the dnality transformation. This snggests 
an idea that if one formnlates the axion theory in terms of the antisymmetric gauge held, then 
the low mass of the axion will be protected not by global bnt by local (i.e. gange) invariance, 
and it will not be destroyed by qnantnm gravity. 


Of conrse, one may immediately argne that this cannot work. Indeed, if the axion charge can 
disappear in the standard formnlation of the theory, its disappearance may have an adeqnate 
description in the theory of an antisymmetric tensor held. However, this argnment has some 
caveats, since in fact these two theories are not completely eqnivalent at the qnantnm level. For 
example, conformal anomaly associated with the psendoscalar held dihers from the conformal 
anomaly in the theory of the antisymmetric tensor held ^0 . 


There is another problem which appears to be mnch more important in the context of onr 
discnssion. Whereas a massless antisymmetric tensor held in a topologically trivial space can be 
converted into a psendoscalar, not all psendoscalars can be replaced by antisymmetric tensors. 
The necessary condition is that the Lagrangian of a psendoscalar held shonld depend only on its 
derivatives. 


Indeed, if one starts with the theory of the antisymmetric tensor gange held, dnality transfor¬ 
mation between the held strength of the antisymmetric tensor held and derivative of 

the psendoscalar held d^a exists and can be nsed to prove the eqnivalence of these two theories 
Dnality transformation is always possible from the 6 ^i,-side to a-side and one ends np with 
the theory of the massless psendoscalar held with derivative conpling only. Vice versa, if one 
starts with the theory of the massless psendoscalar held with derivative conpling only, one can 
nse the dnality transformation and have an eqnivalent 6 ^i.-version of the theory. 
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If the effective action depends on the psendoscalar field a withont derivatives (and this is the 
case if the axion held has a small mass!), no clear information about the relation of this theory 

to the 6„i,-theory was available. To clarify the relation between bnu- versus a-theories one should 

IsP' 

understand how one can describe appearance of the small mass rua ~ terms of the gauge 

theory of the held h^y. Is this ehect possible at all, or is it associated with some kind of gauge 
symmetry breaking? 


We have analysed this question and found that the axion mass generation can be consistently 
described in terms of the antisymmetric tensor held, and that this ehect does not involve any 
gauge symmetry breaking. This ehect is quite interesting in its own terms, independently of 
the initial goal of our investigation. It provides a generalization of the phenomenon studied 
by Polyakov in three-dimensional compact QED, where the massive scalar excitation appears 
in presence of monopoles In our case the mass of the antisymmetric tensor held appears 
because of its interaction with the usual QCD instantons. The deep physical reason why the 
antisymmetric tensor held can acquire mass without breaking gauge invariance is that both the 
electromagnetic held in d = 3 space-time and the antisymmetric tensor held in d = 4 space- 
time have only one degree of freedom. Thus, the condition of transversality for these theories, 
which typically protects excitations from becoming massive, does not apply to physical degrees 
of freedom in these theories. We will describe this ehect in Section 2 of this paper. 


Even though there is no reason to expect that the Peccei-Quinn symmetry is protected by the 
gauge invariance of the antisymmetric tensor held, one may still hope that the symmetry violation 
should be very small for some other reason. For example, it is not so easy for a quantized axion 
charge to be squeezed into a wormhole or a black hole: They should be large enough to eat a 
unit of the axion charge. In the language of Euclidean quantum gravity this translates into the 
question of what is the Euclidean action of the wormhole which could eat a unit of the quantized 
axion charge. If this action S is large enough, then one could expect that the violation of the 
Peccei-Quinn symmetry is strongly suppressed by a factor e~^. Our estimates of the effects 
related to the term giMp{^ + 4)*) indicate that in order to suppress dangerous effects of global 
symmetry violation in the axion theory one should have the wormhole action S > 190. 


Wormholes which could eat the global charge have been Erst discovered by Giddings and 
Strominger in the formulation of the axion theory in terms of the antisymmetric tensor held [Q . 
The simplest of their solutions corresponds to the fixed value of the radial component of the 
scalar held. The psendoscalar representation of the Giddings-Strominger wormhole was obtained 
by Kimyeong Lee |^. In Section 3 of this paper we will re-derive their expression for the 
wormhole action. Our result for the value of the action of the wormhole conhguration without 
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the boundary at the wormhole neck coincides with the result obtained by Giddings, Strominger 
and Lee. However, we point out that if one takes into account boundary terms including the 
contribution of the boundary at the wormhole throat, the action becomes about three times 
smaller. In any case, the action is proportional to which is as large as 10^ for the axion 
theory with /o ~ 10^^ GeV. Thus could suggest that for the axion theory one does not have 
any problem whatsoever since the symmetry violating effects will be suppressed by the factor 
~ which is more than enough to explain why gi < 10“®^. 


Unfortunately, this attitude proves to be too optimistic. As was hrst pointed out in [15 


m 


realistic models of the axion field the radial component f{x) of the axion field on the wormhole 
solutions does not remain equal to /q. Near the wormhole throat this field typically acquires some 
value of the order of Mp S> /q. A detailed investigation of solutions with an account taken of the 
spatial dependence of /(x) was performed by Abbott and Wise ||^ and by Goleman and Lee 


for the case without spontaneous symmetry breaking. The corresponding Euclidean action which 
was found in these papers linearly diverged on extremely large length scales. It was argued that 
despite the action is infinite, wormholes do lead to charge nonconservation and global symmetry 
breaking since the corresponding effects appear on a relatively small scale, where the large scale 
behavior of the wormhole solutions is irrelevant. 


Unfortunately, the most interesting case of the theories with spontaneous symmetry breaking 

It required some additional work to obtain results in a 


was only briefly mentioned in [15, 16 


form in which one could compare them with the expectations expressed in |T^ . Perhaps this was 


the reason why the authors of ref. |T^ did not make any attempt to use the results obtained in 

i, n, |T5|, m. 


In Section 4 of this paper we describe wormhole solutions in several different theories with 
spontaneous symmetry breaking. Whereas in some cases we could obtain important information 
about these solutions by the methods of ref. [^, in general it was necessary to use numerical 
calculations. These calculations were extremely tedious, especially for the axion theories with 
/o ~ 10^^ GeV. The results which we obtained are in a qualitative agreement with the expectations 
of 1]^ [^. We have shown that for a very wide class of potentials the action is finite and to a 
good accuracy is given by a simple expression S' ~ In This means that if the global symmetry 
breaking is suppressed by e~^, this suppression is approximately given by the factor This 
factor is of the order of 10“^ for the axion theory with /q ~ 10^^ GeV, and it is of the order of 
10“^ for the texture theory with ft ~ 10^® GeV. This is clearly insufficient to save the axion and 
the texture theory. 
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In fact, the situation becomes even more complicated if one takes into account that, according 

to [^, the symmetry breaking vertex operators g-n ,——-i" h.c. are suppressed only 

by the part of the Euclidean action S corresponding to integration in a small vicinity of the 
throat of the wormhole. In this case suppression of the global symmetry breaking in the theory 
—+ A|<h^p practically disappears. 


One could expect that this result should be strongly model-dependent. Indeed, we have 
found that near the wormhole throat the held f{x) typically becomes as big as Mp. In certain 
theories, such as the theory of superstrings or supergravity, the effective potential may acquire 
large additional terms, which exponentially grow at large /. In such theories the behavior of the 


wormhole solutions near the throat becomes quite different from the one envisaged in |15, 16 


Here our use of numerical methods was absolutely crucial. The results, however, have not been 
very encouraging: Even if one considers the effective potentials growing at large / as fast as 
exp the resulting Euclidean action remains quite small. Thus, one cannot make the violation 
of the global symmetries small by changing the effective potential of the scalar held in any 
reasonable way. 


Fortunately, during our investigation we have found several ways to hx this problem. First of 
all, when we make the effective potential more and more steep (keeping its minimum at /o), the 
corresponding action tends to increase towards the very large action of the Giddings-Strominger- 
Lee wormhole. We have found that the main reason why it happens is an increase of the size 
of the wormhole: The action is (approximately) proportional to the square of the radius of the 
wormhole throat i?(0), S ~ M|i?^(0). According to our calculations, increasing of the minimal 
radius of the throat just a few times as compared with can make the “natural inflation” 
scenario viable. The situation with axions and textures is more complicated, but still these two 
theories can be saved by an increase of the radius of the throat up to about lOMp ^ or 15Mp^. 

It is very difficult to increase this radius by making the effective potential steep. However, 
there may be some other reasons why the wormhole throat cannot be small. For example, in 
string theory the effective “minimal length” may be somewhat greater than the Planck length 
Mp^. Our investigation contained in Section 5 indicates that the size of the wormhole throat 
can be very large in Kaluza-Klein theories with a sufficiently large radius of compactihcation. 
This suggests that the axion theory can be quite viable in the context of a theory in which 
the gravitational effects on the length scale I < 10 cannot be described by the standard 
Einstein theory of gravity in a four-dimensional space-time. Another example pointing in the 
same direction is related to corrections and effects of conformal anomaly. We show that an 
account taken of conformal anomaly (even if the corresponding terms are relatively small) may 









lead to disappearance of the wormhole solutions. We also discuss the observation made in |^, ^ 
which suggests that there are no wormhole solutions in certain versions of the string theory. 


In Section 6 we discuss a possibility of an additional strong suppression of the wormhole effects 
because of the Gauss-Bonnet term This term does not change any observational 

consequences of the Einstein theory, but it tends to suppress transitions with the change of 
topology. Similar terms appear in the heterotic string theory with the coefficient proportional to 


a'. We show that these terms may suppress the wormhole effects by the factor e W = e ^ 
where g is the gauge coupling constant, is the stringy mass scale. This result is very similar 
to the standard result e for the suppression of the instanton effects in QCD. The possibility 
to have this suppression factor smaller than is quite consistent with the present picture of 
stringy phenomenology. This suppression becomes even stronger if one adds the usual part of the 
action to the topological contribution discussed above. 


Our results and conclusions are summarized in Section 7. A considerable part of our results 
is based on numerical investigation of differential equations for wormholes. In many cases it 


was impossible to solve these equations using standard numerical recipes [^]. In Appendix we 
describe an improved method which we have used in our work. 


2 Axion Theory. Pseudoscalar versus antisymmetric ten¬ 
sor. 

2.1 Pseudoscalar formulation of the axion theory 

The simplest version of the pseudoscalar axion theory []^ adds to the standard model Lagrangian 
CsM the following terms 

£ = + (9 + /■“<*>' . (3) 

For a review of various models see 0. Peccei-Quinn global U{1) symmetry 

a{x) —> a{x) + C (4) 

is broken spontaneously when the operator has a non-vanishing vac¬ 
uum expectation value The term can be represented as a total derivative. 
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= dfj,K^, where 


( 5 ) 


K,{x) = - If^A^AZ) . 

The fact that F^"^^ is a total derivative is not sufficient for providing the symmetry since the 
action is not invariant in presence of instantons. The variation of the action with an account 
taken of the surface terms is 


C 


SSeff = -r 
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fo 327r2 




Cj^ 

fo 327r2 


d^xd^ ^ 0 . 


( 6 ) 


Spontaneous breaking of the global 17(1) Peccei-Quinn symmetry allows to generate a potential 
for the axion held. This gives the axion a small mass ~ The lightness of the 

JQ JQ 

axion is provided by the fact that fo S> Aqcd- 


Note, that the Peccei-Quinn symmetry is global. If the axion would not interact with non- 
Abelian helds, this symmetry could be promoted to a gauge one. Indeed, one can introduce an 
Abelian vector held p and replace everywhere by c?^a-|-eA^. Then the global PQ symmetry 
becomes the local one. 


a —>• a + C{x) , 

^ A^- d^C{x) . (7) 

However, in the presence of a non-Abelian gauge coupling, which is an essential part of all realistic 
axion models, one cannot promote the global PQ symmetry to the local one. Indeed, the non- 
Abelian coupling can be represented in the form d^aK^^. After the promotion of the global PQ 
symmetry to the local one we would obtain the term eA^K^ in the action. This term would violate 
the non-Abelian gauge symmetry since the variation of the Chern-Simons term vanishes only 
when it is coupled to the longitudinal part of A^ but not to the transverse part of it. 

Still it is possible to represent the global PQ symmetry as a local one if one goes to a dual 
formulation of the axion theory, in which the axion pseudoscalar held a is represented by the 
antisymmetric tensor held This version of the axion theory naturally appears in the context 
of string theory 

2.2 Dual version of the axion theory with the non-Abelian coupling 
in the broken symmetry phase 

The non-interacting two-index antisymmetric tensor gauge held hp^ with the gauge symmetric 
action \{d\yhpfj]Y is equivalent to a non-interacting massless pseudoscalar held with the action 
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^{dfj_ay. The simplest way to see this is to perform a duality transformation of the classical action. 
This leads to a replacement of d^a by a pseudovector which is dual to the tree tensor 

held strength d[^bpf^] of the 2-form held bp^j. Another way to see this equivalence is to perform 
gauge hxing of the gauge symmetry 5bp„ = 5[pAo-], SA„ = d„A, which requires two generations 
of ghosts. In addition to 6 fe^j^-helds one gets 8 anticommuting ghosts in hrst generation and 3 
commuting ones in the second generation. This provides a net number of propagating degrees of 
freedom equal to one commuting held. 


The equivalence of the massless pseudoscalar and massless antisymmetric tensor held coupled 
to the non-Abelian gauge held has been shown in starting with the hrst-order type action 
using duality transformation. In what follows we will perform this duality transformation for the 
theory which describes the complex scalar theory with the axion 6 coupled to the non-Abelian 
vector helds, 

c = \dp<!>\^-vm)-dpeQ ^. ( 8 ) 


Here 


^ a: (Ff’ - ‘^h^A'tAD = 


32ir2^"'‘ ' 


(9) 


One can also write this Lagrangian as 

- d,e ni- + £(/). 

where 


and 


£(/) = - V(f) , 

* = exp(^) . M , 


( 10 ) 

( 11 ) 

( 12 ) 


However, it is useful to start with a more general Lagrangian which depends both on the 
pseudoscalar 9 and on some pseudovector held H^, 

Ce,H = ^^p9{H^^ + znn + + £(/) . (13) 

One can solve the held equations for Hp, 

Hp = -tfdp9 , (14) 


and on shell for Hp the Lagrangian acquires the form (pTf ), which describes the pseudoscalar axion 
held 9. 
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( 15 ) 


On the other hand, one can vary Cq^h (0) over 6{x) and obtain the constraint on 

,2 




The solution to this constraint is 


— _ TJ 


Hupa — d[,ybpa] + 327]-2 ~ abc^^pA'a]) i 


(16) 


where 5[,y6po-] is the held strength of the 2-forni held. 


The Lagrangian which follows from (^) describes the antisymmetric massless held hp,^ in¬ 
teracting with the non-Abelian vector helds as well as with the radial component of the scalar 
held, 

(17) 


Pk. = + £(/) ■ 


There are no global symmetries in the dual version of the axion theory, but there are two types 
of local symmetries, the Maxwell-type gauge symmetry of the 2-form held. 


dbpu = , (5A^ = dyk , 

and the Nicolai-Townsend-type non-Abelian gauge symmetry 


5Al = , 6bpy = 2^A^d[pA,y] 


(18) 


(19) 


Will these local symmetries of the dual version allow us to avoid the problem which destroys 
the nice features of the axion theory in the standard version? To address this issue we will hrst 
study the origin of the light axion mass generated by QCD instantons in the dual version of the 
theory.^ 


2.3 How does the axion become massive in the dual theory? 


The possibility that a scalar or pseudoscalar particle can acquire a non-vanishing mass due to 
non-perturbative ehects does not look very surprising. However, in the dual formulation of the 
axion theory the axion is massless because of a gauge symmetry ([l^). It is commonly believed 


'^We are grateful to M. Dine for the suggestion to investigate this problem. 
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that gauge symmetry protects massless particles from becoming massive. This would solve all 
our problems. However, if it were true, then the axion field in its dual formulation would not get 
its mass nia ~ Therefore before going any further we must hrst resolve this puzzle. 

With an account of QCD instantons the effective action of the pseudoscalar field acquires a 
potential, 

^a = ^{d^,aY^ + A^QCDCOs^. ( 20 ) 

Starting with this action one cannot make a transition to the dual version of the theory with the 
-field using any of the methods discussed above. For example, one cannot perform duality 
transformation of the action starting with the action depending both on and as we did 
before. Indeed, there is now a term in the pseudoscalar version of the theory depending on a 
rather than on d^a and the procedure does not work. If we will try to go from the side of the 
gauge invariant 6 ^ 1 / -theory and perform duality transformation, we will get no terms depending 
on a, but only terms depending on cl^a. 

Does the impossibility to perform duality transformation mean that the dual version of the 
axion theory is incapable of explaining the mechanism of generating a small axion mass? Will 
the same mechanism protect axion from getting very heavy? 


To understand the effect at the level which is more subtle than just performing duality trans¬ 
formation over the effective action we will hrst go back to the Polyakov model of compact QED 


in d = 3 |22|. In this theory the massless scalar is dual to the massless Abelian vector and 


not to the -held. 


Polyakov’s main idea is that the existence of magnetic monopoles, which are instantons in 
d = 3, changes the behavior of the correlators: due to instantons there are no massless transverse 
excitations in the system. Instead there is an excitation corresponding to a longitudinal compo¬ 
nent of a gauge invariant operator, which may be interpreted as a massive scalar. The ehect has 
a non-perturbative nature and does not violate gauge symmetry. 


We will show that something very similar happens with the antisymmetric tensor held in our 
d = 4 theory. Let us hrst describe the situation in d = 3 [^. One starts with the dual to the 
vector held strength. 


Ha 


( 21 ) 


If one allows certain magnetic charge density in the system of the form p{x) = X) — Xa), in 
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the quasiclassical approximation this would correspond to 


2mk^ 


( 22 ) 


The main steps in the calculation of the correlator of two gauge invariant operators are the 
following. The bare part without monopoles is 

- '^) . (23) 

This expression corresponds to the free photon propagator coming from the action 
addition, there is a contribution from the instantons such that the full correlator is 

{H,(k)H,(-k)) = {H,(k)H,{-k)Yk) + {p(k)p(-k)) . (24) 


The correlator of charge densities is calculated from the generating functional for the charge 
density of the plasma and is given in [^, 


Z{rf) = (e* ^ 


Z(0) 


Vx exp - 


27r/ 


(V(x - ri)f - cosx 


(25) 


where exp , C = 0(1) is some constant. The second variational derivative over r] 

gives the expression for the correlator of p, 


{p{k)p{-k)) = 


-V 

27rJ 


k^ 


M^ + k'^ 


/e\2 

V^j M‘^ + k^ ■ 


(26) 


The hnal answer is 


:{H^{k)H,{-k)) = 6^, - 


k^ky 

M^ + k‘^ 


(27) 


Thus the correlation function of two gauge invariant operators has only one longitu¬ 

dinal excitation as if one would calculate the two-point correlator of the derivatives of a massive 
scalar held. 


For the antisymmetric tensor held in d = 4 things work not exactly the same way but very 
close to it. The only gauge invariant operator in our theory where we keep the antisymmetric 
tensor coupled to the non-Abelian vector held is 


ppp — _ 

2 




9 


327r2 


a:{f, 


apcr 


|/al,o4lj2lJ) 


= . 


(28) 


This is the only operator which is gauge invariant under both Maxwell and Nicolai-Townsend 
non-Abelian symmetry. Separately, the operator which is dual to a held strength of the 
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antisymmetric tensor field is invariant nnder the Abelian gange symmetry, bnt it is not invariant 
nnder the Yang-Mills symmetry withont the Chern-Simons term. 

To examine onr problem we may ignore small qnantnm flnctnations of the radial component 
of the scalar held. These hnctnations are irrelevant when one investigates the possibility that 
gange symmetries protect axion theory from getting a mass. Thns, the action we consider is the 
action ( 0 ) with the hxed radial component, /(x) = /o, 

• ( 29 ) 

We may start by treating the conpling of the to the non-Abelian held pertnrbatively, 

i.e. we may stndy the correlator of two gange invariant operators by treating the conpling 
as a correction to the valne of the correlator of two operators withont conpling. The bare 
correlation fnnction of two operators mnst be transverse since 

= 0 . (30) 

Thns, as in the Polyakov case, we may expect 

- ^) + '^(p[k)p{-k)) + ... , (31) 

where the dots correspond to the corrections to a transverse part of the correlator. We have 
parametrized the longitudinal excitation by some 2-point correlator of a “charge density” p{x). 
The divergence of our gauge invariant operator which we may associate with p, is given by 

fpd,H'‘(x) = = -ipix) . (32) 

The correlator of two divergences of the Yang-Mills Chern-Simons currents was calculated by 
Shifman, Vainstein and Zakharov in ref. |^. They have shown that in a certain approximation 
it is proportional to the axion mass 

F“](-*:)) ~ ml . (33) 

This serves as an indication that one may obtain in the 6^jy-theory the longitudinal excitation 
of the gauge invariant operator corresponding to a scalar massive particle instead of a massless 
bfxu -held. 


One can conhrm these expectations by calculating the correlator {H^{k)Hj,{—k)) directly, 
without separating our operator into its free part and the interaction part For this 
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purpose we will perform the generalized duality transformation in the functional integral describ¬ 
ing the theory. Consider the “hrst order” functional integral for the action we have considered 
above with the additional term including the source r]^ to the current 


Z{rf) = J VaVH^VA'^ exp 


* / + Cym) 


(34) 


The functional integration is performed over the pseudoscalar a as well as over the pseudovector 
and over the Yang-Mills helds. For the sake of simplicity we do not write down explicitly the 
integration over the ghosts related to the gauge hxing of the non-Abelian symmetry. We may 
evaluate this functional integral by integrating over the pseudoscalar a hrst. This will produce a 


constraint given in eq. (IT). In this way we will get rid of the a-integration and the remaining 
integral will reduce to the integration over the constrained Equivalently the integration over 
the constrained may be replaced by the integration over the unconstrained 6 ^ 1 . together with 
the proper gauge-hxing procedure in the path integral. As before, H^{b, A) is dehned in terms of 


and Yang-Mills helds in eq. (PT]). Thus we get 


= /' 


= / VbyyDCYt.'DAl exp 


' / A) + Cym + Aax + •Cja 


(35) 

where we have written down the integration over the complete set of ghost helds Cgh related to 
both types of gauge symmetries. £gfix is the gauge-hxing part of the action and £gh is the action 
of the ghost helds. By varying this functional twice over we will get the correlations function of 
two gauge invariant operators A) dehned in eq. (pT[) and calculated in the dual version 

of the axion theory. 


On the other hand, we can perform the integration over hrst. This does not change the 
fact that the second derivative of Z{j]) gives the correlator of two It is just an alternative 

method of calculations. After integration over we get 


Z{r]) = / VaVAl exp 


--{id^a + ri^y + 


g‘^a 


327r2/o ^ 


(36) 


The integration over the instanton conhgurations in the dilute gas approximation produces the 
ehective action for the held a 


Ari) 


VaVAl exp 


/ I Qj 

-(8^a - iri^y + AJcd cos -y + Cym 


(37) 


Double variation of this generating functional gives us an expression for the correlation function 
of two operators in the following form: 


U\H,{k)Hy{-k)) = 5 


fll/ 


kf^ky 

m? + 


(38) 
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This result is in a complete correspondence with the Polyakov result for d = 3. The hrst term in 
the r.h.s. of eq. (|38D comes from the term in the generating functional, and the second term 
comes effectively from the correlator of two derivatives of the massive field a. 

Our conclusion is the following. The antisymmetric gauge held is coupled to the Yang-Mills 
helds. With an account taken of instantons the behavior of the system becomes that of the 
massive scalar held theory. Note that this is not a mass of the held Rather it is the mass of 
the only gauge-invariant degree of freedom associated with 

Thus we have found a mechanism of getting small QCD mass for the axion in the dual version 
of the theory, which does not violate any of the gauge symmetries. This result has an important 
implication that there is no reason to expect that the gauge symmetry which exists in the dual 
formulation of the axion theory can protect the axion mass from becoming very large because of 
the gravitational effects. 


3 Axionic instantons, or wormhole solutions with fixed 
radial component of the scalar field 


The axionic wormholes which may provide mass to the axion via gravitational effects have been 
found originally in the dual version of the axion theory with antisymmetric tensor field 0. The 
theory of the complex scalar held does not have such wormhole solutions unless the functional 
integral is supplemented by the proper boundary conditions. The corresponding investigation has 
been performed in ||^ and in [|1^. The conclusion was that there exists a consistent procedure 
to obtain the same wormhole solutions in both versions of the axion theory. 


In this section we discuss wormholes with a frozen radial component of the held, like Giddings 
and Strominger ^ and Lee [^. We will mainly reproduce their results. However, in addition we 
will discuss the subtlety related to the boundary terms and the value of the Euclidean action. We 
will hnd out that if one considers the conhguration with the boundary at the wormhole throat, see 
Fig. 1 (a), and calculates the contribution to the action from the boundary terms in a standard 
way, one ends up with the action which is only (1 — |) ~ 0.36 of the original action, which was 
23| without an account taken of the boundary terms on the wormhole throat .0 


calculated in 


^We are grateful to A. Strominger for the discussion of this subtlety. 
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Giddings and Strominger have found their wormhole solution in the string-inspired theory 
with the action 

M ,^S(K - A-,) . (39) 

Here = d\^^hyx] is the held strength of the antisymmetric tensor held. The last term in the 
action is the Gibbons-Hawking surface term. 

Meanwhile, Lee obtained wormhole solutions in the theory of a pseudoscalar axion held 6 = -^ 
with the action 

®- § L ■ <«) 

The wormhole geometry has the following form 

ds^ = dP + R{rYd^Vt^ . (41) 


As we already mentioned, there are no true solutions of the Lagrange equations following from 


(^). However, it was pointed out in that these solutions appear if one takes into account 

the charge conservation condition in space (|4l|) . The global charge n, dehned as an integral over 
the 3-space from the zero component of the Noether current. 


n = y Wp9' = f^O' = const , 


(42) 


IS an 


integer upon quantization.]^ 


The conclusion of ref. ||^, was that if one makes the variation of the action under the 
condition n = const, one obtains equations of motion which are equivalent to the equations in 
the theory ([39|). In particular, the gravitational equation of motion in both theories looks like 




(43) 


Here i?(0) is the size of the throat of the wormhole dehned by the condition = 0. It is given 
by 


fi(0) = 


n 




(44) 


^Note that in the situation with spontaneous symmetry the global charge is carried not by charged particles, 
but by the vortices of the classical scalar field with the time-dependent phase 9 . 
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The wormhole geometry R{r) can be expressed analytically in terms of elliptic integrals 0. 


The wormhole action with account taken of the boundary terms both at the outer boundary 
at r = To, To —»• cxD, and at the inner boundary at r = 0 is 

S A/f^ c 

5'totai = J dvt drR{r)R'{r){l - R'{r)) . (45) 

The integration can be performed as follows 


‘S'total 


3M| 

dvr 


271^ - 


2 j_r2(o) 


dR^ 




(46) 


Giddings-Strominger-Lee on shell action (without the inner boundary terms) is 




nobound 


37rM| 


37r^ 


drR{r){l - [i?'(r)]2) = —M^R^{0) 


The boundary term which we added, 

^ho„nH = r drR{r){l - R'{r)) = 

JO 


377^ M^R\0) 


(47) 


(48) 


2 Ja 

is —0.637 of the action without the boundary term. To understand better the boundary term 
contribution consider it in the form 

‘S'bound 


3M2 


Svr 


dQ / dr 


A 

dr 


R^{r){l- R'{r)) 


(49) 


The contribution comes only from the throat where R{r) = R{0), R'{r) = 0. It equals 


5, 


bound = -^2w^R^(0) = -^M^fl=(0) , 

OTT 4 


(50) 


which is in a complete agreement with ([481) . Note that the extrinsic curvature term K at the 
boundary at r = 0 (the second term in eq. (|49|) ) vanishes, in agreement with the discussion in 
p6| , but an additional non-vanishing contribution comes from the term Kq, which was added to 


remove the divergence of the action at the outer boundary. 


The total action for the configuration with the boundary at the throat is thus only 0.363 of 
the value of the action without the surface term. This gives an additional support to the idea 
that in gravitational problems one has to be very careful with the boundary terms. The outer 
boundary surface term K has to be corrected by Kq, otherwise the action is inhnite. Does it 
mean that the surface term has to have the same functional form K — Kq on both boundaries? 
If the answer is yes, we have to subtract ^ part of the action and make it almost three times 
smaller than the action obtained in P, 
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One should note, that the problem of obtaining a proper contribution from the inner boundary 
is very complicated. For example, recently it was argued that an adequate account taken of the 
inner boundary of an extreme charged black hole changes its Euler number |^, which makes its 

However, even for the black hole case this issue is rather nontrivial. 


entropy zero 
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The situation with the wormholes is even more complicated and ambiguous. If one considers 
the wormholes without the inner boundary contribution as in P, |^, |T^, one has to add 

to our result for the action the universal term ^Mpi?^(0), which depends only on the size of 
the wormhole throat. One may also want to calculate the action on a symmetric conhguration 
—cx) < r < + 00 , Fig. 1 (b). In this case one will not have any contribution from the inner 
boundary, and the action will be two times greater than the GSL action 


In what follows we will work with the action dehned with both boundaries and with the surface 
term K — Kq on both boundaries since this prescription seems to be more internally consistent. 
Another advantage of this prescription is that it gives us the smallest action as compared with 
other prescriptions mentioned above. Therefore if we hnd a way to avoid the strong violation of 
global symmetries in our approach, we will simultaneously solve the corresponding problem in 
other approaches as well. 


The radius of the throat of the Giddings-Strominger-Lee wormhole depends on the parameter 
/o and on the value of the charge n. We may therefore express the action in terms of these 
parameters as follows: 






(1 


2 \ nMp 


TT 


/o 


For completeness, we will give here also the value which has been obtained in 
account taken of the inner boundary. 


(51) 
231 without an 


5 , 


nobound 



nMp\ 

~) 


(52) 


If one takes the smallest action (^) with n = 1 and /o = 10^^ GeV, one obtains an enormously 
strong suppression ~ exp(—10®). This would immediately solve the problem of the global sym¬ 
metry violation. Unfortunately, however, things are much more complicated. As it was pointed 
out in |T3|, it is almost impossible to keep the held f{x) close to /o on the wormhole solution. In 
what follows we will show that if one allows the held / to depend on r, this held typically grows 
to / Mp 3> /o near the wormhole throat. Therefore one needs to make a separate investigation 
to calculate the wormhole action in realistic theories with spontaneous symmetry breaking. This 
investigation will be contained in the next section. 


20 























4 Wormhole solutions with dynamical complex scalar 
field and spontaneous symmetry breaking 

4.1 Equations for the scalar field in the wormhole geometry 

In this section we will discuss wormhole solutions for several different theories with spontaneous 
symmetry breaking. In particular, we will investigate the dependence of the action on the value of 
the vacuum expectation value /o of the radial part of the scalar held far away from the wormhole. 

We will study interaction of gravity with the complex scalar held <I>(a:) = The 

corresponding action is 

S = I + + V{f)^ - ^ Ko^) . (63) 

We will assume that the potential V{f) has a minimum at some value / = /q. The vacuum 
energy vanishes in this minimum, 

^'(/)l /=/0 = n/)l/=/o = 0. (54) 


Equations of motion corresponding to the analytic continuation of the Euclidean theory of 
the 2-form version of the axion theory |^, |^, are 

3Rf dVif) 


r 

Stt 


R df 

2 f>2^ 


n 




= 0 , 


= 0 


(55) 

(56) 


and the value of 6' = been already substituted in the equations. 

Using these equations, one can derive the following expression for the wormhole action: 


27r 


dr + 


- R) 

dvr 


(57) 
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We are looking for the wormhole solutions with the geometry given in eq. (^T]) and with the 
helds R{r) and /(r) solving the system of equations (|^), (|^). Our boundary conditions are 


R'{0) = 0 , (58) 

/'(O) = 0, (59) 

/'(r) ^ 0 , r ^ cx) , (60) 

f{r) /o , r ^ oo . (61) 

We will examine the following potentials: 

Vdf) = jif-fR. ( 62 ) 

n(/) = ^ (/'- 3/„V" + 2 /„') . (63) 

n(/) = (64) 


The hrst potential is a standard potential of a theory with spontaneous symmetry breaking, the 
second one is inspired by some phenomenological models based on supergravity with /q ~ ^ 

the third one is inspired by string theory. We will hnd the wormhole solutions numerically using 
the fact that the held / has a dehnite value /o at inhnity. In particular, for axions we will 
be interested in /q ~ 10^^ GeV and for the textures in /q ~ 10^® GeV.^ Typically the radial 
component of the axion held / will grow many orders of magnitude from the value /o to / ~ Mp 
when approaching the wormhole throat. 


Note that our expression for the total action ([57|) does not have any explicit dependence on 
the ehective potential V{f) and on the charge n. It is important also that 0 < R'{r) < 1 for all 
wormhole solutions which we are going to study. As a result, the integrand is always positive, 
despite the fact that the gravitational contribution to the action may be negative. One of the 
consequences of this result is that if we make a cut-oh at some radius Tc and integrate from r = 0 
to r = Tc, the result will be always smaller than the total action 


4.2 Wormholes in the theories without symmetry breaking 

In this paper we study the theories with spontaneous symmetry breaking. However, on a small 

scale R < the wormholes in the theories with the effective potential 

"'^We should emphasize again that the value of the parameter /o may be quite different from 10^^ GeV. However, 
as we will see, our results depend on /o only logarithmically. Also, for textures one should take a theory with 
another group of symmetries, which, however, should not considerably change our results. 
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described in ||T5| , are very similar to the wormholes which we have found in many theories with 
spontaneous symmetry breaking. Therefore we will briefly discuss here the wormhole solutions 
in the theory following pr5[ . 


Near the wormhole throat the solution was obtained in numerically. This solution should 
be matched to the solution obtained analytically at large r. According to eq. (|5^ , far away from 
the throat one has R{r) = r. Therefore equation (|^ in this case (for n = 1) reads 



1 

4^4y3^6 


0 . 


(65) 


This equation has an exact solution 



where 6 is determined from the equation 


( 66 ) 


1 

47r^d^ 


-AS^ 


1 . 


(67) 


If one adds to the effective potential the term with > 0 (no spontaneous symmetry 

breaking), the asymptotic behavior at large r changes to / ~ (27r^mr^)“^/^. 


It was assumed in that the contribution to the action from the region near the wormhole 
throat (until the solution approaches its regime (|66l)) is very small. On the other hand, integration 
far away from the throat gives the contribution 


S = 2tt^ / r^dr 


l /'2 ^^- 

2'^ 


+ VU) 


( 68 ) 


The total action integrated up to some rmax has the following general form: 


3A 

5total(j’max) ^ “(1 “ ^) ln(mr„) + mr^ax + A5 . 


(69) 


Here corresponds to the place where the numerical solution near the wormhole throat matches 
the solution (|66|), AS" stands for several other terms which have not been determined in |]^ 


One can calculate AS using our expression for the action (|F7|); typically this term is fairly small, 
A^ = 0(1). 


Thus the two most interesting terms in (|6^) is the logarithmic term and the term which linearly 
diverges at large rmax- As we will see, a similar logarithmic contribution appears in the theories 
with spontaneous symmetry breaking as well. However, the linear divergence at r^ax —>■ cxo is 
a particular property of the theory without symmetry breaking, which is a consequence of the 
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asymptotic behavior / ~ at large r. One could conclude that the wormhole action 

is inhnite in the theories without spontaneous symmetry breaking, and therefore these theories 
cannot lead to the global symmetry violation. 


However, it was argued in [|^ that this is not the case, and in fact the effects of global 
symmetry violation are quite signihcant. It was suggested that these effects are suppressed only 
by some small part of action Syj coming from the region close to the throat of the wormhole. 
The value of was not calculated in Q, but it was estimated to be 0(1). We will return 
to the discussion of this issue when we will consider wormholes in the theories with spontaneous 
symmetry breaking. 


2\2 


4.3 Wormholes in the theory with the simplest potential — /q) 


Equation for the scalar held / in the theory |(/^ — /q)^ on the wormhole conhguration looks 


very similar to eq. (|6^: 


^ - Wf - fS) - 


= 0 . 


(70) 


Therefore the wormhole solutions at / 3> /o behave just as their counterparts in the theory 
without spontaneous symmetry breaking. In particular, far away from the wormhole throat 
/ = ^, where 6 is determined by ^^ 1^4 — = 1. 


The main difference between the wormhole solutions with and without symmetry breaking 
appears on the scale r > ^4A/2vr^A^ /q”^, where the held / approaches /o very rapidly (though 

not exponentially, as anticipated in |^). We have found that at large r 


f{r) - fo = 


Stt^A/qV 


5.^6 ’ 


(71) 

-V5 

0 


which leads to a hnite (and very small) contribution to action from the region r > (4:^/2%^f^ 
Thus, in our theory we do not have any problems with inhnite wormhole action. 


To study the behavior of our solutions in all regions from r = 0 to r —> cx), and to avoid having 
numerical uncertainties associated with the value of action coming from each of the regions we 
performed a numerical investigation of the wormhole solutions in our theory. This investigation 
was rather complicated since the solutions happen to be extremely sensitive to the boundary 
conditions. Whereas it was possible to obtain some results for /q ~ 10^® GeV using standard 


numerical recipes it was necessary to substantially improve these methods in order to study 


24 
















the physically interesting regime with /o ~ 10^^ GeV. Since this improvement may be useful not 
only for hnding the wormhole solutions, we will describe our method in the Appendix. 


We have found it useful to make our numerical calculations in dimensionless variables p, A(p), 
F and U{F), where 


p = rMp W — , A = RMp \ — , 
V OTT V OTT 


^"MpVT’ t'W = VVW . (72) 


We also introduce the combination Q = 


Equations of motion in these variables are given by 

3A'{p)F\p) dU{F) 2Q^ 


F"{p) + 


^(P) 


A-(p)-l + A(p)M[/(F) + 


dF F3(p)yl6(p) 

_ 1 ..22 A 

F\p)A^{p) 2 


= 0 , 


= 0 . 


(73) 

(74) 


The on shell action is 


poo p 

S„., = ^l dp [A^(p)F'‘(p) + 2A(p)A'(p){l - A(p) 


(76) 


Our numerical solutions of the system of differential eqs. (|7^) for the theory — /q)^ 
depend on the value of the coupling constant A and on the asymptotic value of the held /q. In 
all cases we are interested in the strongest possible violation of global symmetries by gravity and 
therefore we will consider the smallest value of the charge n = 1. We will present most of our 
results for A = 0.1 and A = 1, but in all hgures we will show only the case A = 0.1. Different 
asymptotic values of the held /o are considered, from /o = 10^^ GeV to /o = 10 ^^ GeV. The 
corresponding values of the dimensionless held Tq are 


Fo = 



fo X 2.4 X 10“^^ GeV^ . 


(76) 


For example for /o = 10^^ GeV one has Fq = 2.4 x 10 and for /q = 10^® GeV one has 
Fo = 2.4 X 10 3. The solution for the dimensionless radial component of the scalar held F is 
represented on Fig. 2 where the value of logj^gF is plotted as a function of log^gp for diherent 
seven values of Fq, corresponding to /o = (10^®, 10^^, 10^®, 10 ^®^ 10 ^®^ 10^^) GeV. One can 

see from Fig. 2 that all solutions with diherent asymptotic values of /o behave in the same way 
near to the wormhole throat. This means that the axion type held F which far away from the 
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wormhole was ~ 10“^, increases seven orders of magnitude near the throat to become ~ 1, and 
the texture-type held F which far away from the wormhole was ~ 10“^, increases three orders 
of magnitude near the throat to reach the same value ~ 1 corresponding to / ~ Mp. Thus, the 
solutions which we obtained differ very much from the Giddings-Strominger-Lee solutions with a 
hxed value /(r) = /q. On the other hand, the fact that the solutions at small r do not depend on 
/o conhrm our expectations that the behavior of these solutions near the wormhole throat does 
not depend on spontaneous symmetry breaking. 

This conclusion becomes even more obvious if one consider Fig. 3, which gives the value of 
the function A{p) (i.e. R{r)). This function, describing the geometry is completely insensitive to 
the asymptotic value of the held Fq (i.e. of /o). Therefore we have only one curve for all seven 
cases above. 


One can also express our results in the usual dimensional variables. One can show, in partic¬ 
ular, that if the coupling constant A is very small, then, just as in the theory without symmetry 
breaking [^, the value of the held / at the wormhole throat and the radius of the throat are 
given by the simple expressions independent on A and /q: 


/(O) « Mp 



R(0) « . 

Y3\/27t 


P7) 


Meanwhile the total action does depend on /o, and this dependence is pretty simple. Fig. 
4 represents the value of the action as the function of — InFo for A = 0.1. The black line and 
dots give the total action with the boundary term, the grey one shows the action without the 
contribution of the boundary term at the inner boundary (see the discussion of this possibility 
in the previous section). All data £t the following simple equation for the action 

5'totai = a - & In Fo . (78) 

The values of a, 6 for two values of A are 

a = 0.186 ±0.005, 5 = 1.001 ± 0.005 , A = 0.1 , (79) 


a = 0.188 ±0.01, 5 = 1.010 ±0.002 , A = l. (80) 

Thus the dependence on A is not strong and the total action is small. 

If we consider the configuration without the boundary term, the action increases slightly. The 
values of a and b for two values of A are 

a = 0.850 ±0.005, 5 = 1.001 ± 0.005 A = 0.1 , (81) 
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1.33 ±0.01, 


6 = 1.010 ± 0.002 A = l. 


a = 


(82) 


However, even in this case the total action remains of the order of 15 for /o 10^2 GeV. This is 
much smaller than what we need. 


This conclusion may seem somewhat unexpected. Indeed, the only natural length scale in 
the theory of gravity is Mp How could it be possible for a wormhole with a throat of a radius 
~ Mp ^ to eat a vortex with f^9 7 ^ 0 or a particle of a typical size m~^ 3> Mp Indeed, if 
wormholes do not change the value of the scalar held /, such processes are extremely strongly 
suppressed, as we have seen for the case of the Giddings-Strominger-Lee wormhole. However, in 
our case the total action is rather small, and it depends on /o only logarithmically. A possible 
interpretation of our results is the following. It does not cost the wormhole almost any action 
to squeeze the vortex to the size r A Later (in the Euclidean time r) by increasing 

the scalar held / the wormhole squeezes the vortex to the Planckian size and easily swallows it. 
One may say that our wormholes have a small throat but a very big mouth; they compactify the 
charge before eating it. 


Note that the total action provides the maximum value of suppression of the violation of 
global symmetries. However, as we have already mentioned, this suppression in fact may be even 
much weaker, if it is controlled not by the total action but only by the contribution to the action 


from the small vicinity of the wormhole throat |]^ |I| 


Indeed, nonperturbative effects are controlled by the total action only if one can use the dilute 
gas approximation and consider contribution of each wormhole separately. If the wormholes are 
very compact and their action is very large, then it is indeed the case. Otherwise one may consider 
a possibility that the wormholes carrying away opposite charges can screen the large-scale “tails” 
of each other, and their effective action then will be determined by the integration from r = 0 to 
Tc, where 2rc is a typical distance between the throats of different wormholes. If is not much 
different from the radius of the wormhole throat i?( 0 ), the corresponding action should be small, 
and there will be no suppression of the wormhole-induced global symmetry violation. 


To make an estimate of the cut-off radius (assuming that we are already in the regime 
R ^ r) one may write an approximate condition implying that there are no wormholes within 
the distance 2 rc from each other: 
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Here e appears due to the exponential suppression of the wormhole-like fluctuations on 

the scale Tc, and is our estimate for the subexponential factor. 


Using these results and eq. ( |55D one can obtain the value of which should be 

associated with the effective coupling constant of the operators violating global symmetries. In 
the case we are considering right now it is a pretty easy problem to solve. Indeed, at r ^ -R(O) 
our wormhole solution enters the regime R = r, f = and its action at this stage with a very 
good accuracy is equal to In(compare with (^)). Therefore our condition (|83[) in this case 


reads 
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It is clear that this condition cannot be satished for r > i?(0). Thus, Vc < .R(O) for the wormholes 
in the theory |(/^ — /q)^. This gives S'totai(’"c) = 0(1), which implies existence of the vertex 
operators gMp(^ -|- $*) violating global U(l) symmetry with the unacceptably large effective 
coupling constant g ~ = 0(1). 


Thus, spontaneous symmetry breaking per se does not imply any suppression of the wormhole- 
induced violation of the global symmetries. One should be warned, however, that this conclusion 
was based on a rather crude estimate of using eq. (|^). This equation is based on 

the assumption that if one has many wormholes at a distance 2rc from each other, the action of 
each wormhole is (approximately) equal to the action of a single wormhole solution integrated 
from r = 0 up to Vc- Note, however, that the scalar held f aX r = Vc remains extremely large, 
/(rc) ~ Mp S> /o- Thus one could argue that if space were populated by many wormholes 
displaced at a very small distance from each other, this would not describe our original situation 
where the average amplitude of the radial component is equal to /o -C Mp. This argument does 
not really invalidate the multi-wormhole scenario. The radial part of the held may be quite 
large in the vicinity of each wormhole, but the presence of the charge implies that the phase 6 
depends on r and can be diherent for each of the wormholes. Therefore, even though the value 
of the held near each of the wormholes is of the order of Mp, the average held $ in the whole 
space (after taking the average over huctuations of 6) can have a very small radial component 
/o -C Mp. 


To return to a more solid ground, one should note that even if one does not want to consider 
this multi-wormhole picture and decrease the wormhole action by introducing a cut-oh at r = Vc, 
one still has a problem. Indeed, in the theory with spontaneous symmetry breaking, unlike in 
the theory without symmetry breaking considered in |T^, M, the total wormhole action is hnite 


and small. The largest action which we have obtained for the axion theory with /q ~ 10^^ GeV 
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is only about 15. This is more than ten times smaller than the action S ~ 190 which is necessary 
to save the axion theory. Therefore we will study now other, more complicated models, where 
one may hope to obtain larger values of S. 


Note that if we hnd the wormholes which have a very large action S > 190 given by the 
integration in a small vicinity of the wormhole throat (and we will hnd such solutions), then our 
equation (|8^ will suggest that the cut-off radius Vc should be many orders of magnitude greater 
than i?(0). Therefore in all situations where we will hnd a solution to the problem of strong 
violation of the global symmetries, the ehective coupling constants in our vertex operators of 
the type of giMp{^ -|- $*) will be determined by the total wormhole action Stotai rather than 
by its small part originated by the integration near the wormhole throat. This will eliminate 
all uncertainties with the interpretation of the vertex operators and multi-wormhole solutions 
described above. 


4.4 A more complicated polynomial potential 

The next theory on our list has an ehective potential 

V2(f) = T - 3/o/= + Vt) , (85) 

where /q This potential was suggested to us by M. Dine as a useful phenomenological 

potential for the axion theory which might follow from supergravity. We have found the wormhole 
solution in this theory and performed the calculation of the action. There was practically no 
diherence in the value of the action as compared with the action in the theory — /q)^. 
Again we have recovered the logarithmic dependence on the value of /q. All hgures practically 
coincide with the ones obtained in the previous case. The conclusion is that the mild change 
of the potential responsible for spontaneous symmetry breaking does not change the value of 
the wormhole action and therefore in such theories we have to face a strong violation of global 
symmetries. 

The reason can be easily understood. In the theory — /g)^ the held / far away from the 
wormhole throat behaved as / = ^, with 6 being determined by ^^ 4^4 — A(5^ = 1. When the held 
/ decreases as / = ^, the contribution of the term A/^ to the equation of motion 0 decreases 
as fast as the contribution of all other terms (for / 3> /o). If now one has an ehective potential 
which depends on / as /®, then its contribution to the equations of motion in the regime / ~ ^ 
decreases even faster than other terms, and the held behaves as in the theory |(/^ — /g)^ in the 
small A limit. 
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Thus the existence of the regime / ~ ^ is a very general property of our wormhole solutions. 
This leads to the familiar logarithmic dependence of the action on /q. Therefore it is very difficult 
to increase the action S by changing /q. However, as we will see now, one can considerably increase 
the action if one succeeds to increase the radius of the wormhole throat i?(0). 

4.5 Exponential potential 

The third class of potentials includes an exponential dependence on the radial component of the 
held which forces the held to remain close to its minimum value and not grow so fast as in the 
previous cases, V 3 (/) = (/^ — /q)^. In terms of our dimensionless variables the potential 

is 

UiiF) = (F" - FSf . (86) 

Our previous calculations have shown that the wormhole action is small when the theory allows 
the radial component of the held to grow near the throat. By introducing the exponent to the 
potential we were hoping to achieve several diherent purposes: 

i) We were trying to keep interactions far away from the wormhole the same as in the usnal 
theory of spontaneously broken symmetry. Indeed at small F this potential coincides with the 
standard potential 

ii) When approaching small distances the value of the exponent will change strongly; the 
growth of the held F will not be able to proceed. This should increase the size of the throat of 
the wormhole. 

hi) The exponent in the potential should imitate a gravitational theory which does not allow 
distances smaller than some valnes. In this way we could have a model where the size of the 
wormhole throat cannot be small. As we will see, this will give us large wormholes with the 
action approaching the limiting case of the Giddings-Strominger-Lee wormhole with the radial 
component of the scalar held frozen to its value at inhnity, F = Fq. 

Our expectations were indeed conhrmed by the calculations. We present the wormhole soln- 
tions for hve different values of the exponent (3 = (0,10,100, 300, 500). We work with A = 10“^ 
and /o = 10^® GeV (i.e. Fq = 2.4 x 10“^). Fig. 5 shows the plot of log^QF as the function of 
logiop. inhnity all solutions are given by the same line, corresponding to Fq = 2.4 x 10“^. 
However they behave diherently when approaching the throat. The solutions with the largest 
values of the exponent do not grow as fast as the solution with smaller exponents, the grows in- 
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creases when the exponent decreases. The behavior of the radial component of the held shows the 
tendency to approach the GSL solution where the held F is hxed. In our case this is a property 
of the solution of a system of equations ([73D, ([7^1), whereas the GSL solution is a solution of the 
gravitational eq. (^) only, in which one assumes that F' = 0. It is gratifying that the dynamical 
system with the standard mechanism of spontaneous symmetry breaking can be brought to the 
regime of the almost frozen radial component of the hied which results in the increase of the 
wormhole action and suppression of the violation of global symmetries. The geometric properties 
of the wormhole are represented in Fig. 6, which shows the A{r) and in Fig. 7, which shows 
A'{p) = R'{r). One can see a dramatic change of the geometry with an increase of the exponent. 
With zero exponent we have the same situation as in the theory ^{F^ — Fg)^. The wormhole 
throat is very small, A' becomes close to the flat space value A' = 1 very fast. With the growth 
of (3 the picture changes: A' remains smaller than 1 even far away from p = 0, and the size of 
the throat A(0) becomes much larger. 

The total action was calculated for several different values of the exponent jS. The results are 
plotted in Fig. 8. We have found that in a wide range of values of (3 from 0 to 500 the action is 
given by a very simple expression 

-S'totai ^ a + h(3 . (87) 

The values of a, b for the action with the surface terms are a = 5.7, b = 0.034. The corresponding 
dots and the line in Fig. 9 are black. For completeness we will also give the results for the action 
without the boundary term at the throat (see the discussion of this issue in Sec. H). In this case 
a = 6.8, b = 0.083, see the grey line in Fig. 8. 

Note that the value of the action without the boundary term for the largest value of our 
exponent is twice as large as the one with the boundary term, although with vanishing exponent 
the difference is very small. The reason is that at higher value of the exponent the value of the 
wormhole neck is getting much larger and therefore the term ^M|i?^(0) is also getting larger. 

We have also calculated the part of the action coming only from the area close to the wormhole 
throat, where the geometry is much different from the geometry of a flat space. This part does 
not take into account the logarithmic contribution which comes from the region with R' ^ 1. We 
have called the part of the action which comes from the region from r = 0 to the radius r 
at which R'{r) grows up to 0.9. On dimensional grounds one could expect that should be a 
quadratic function of the size of the throat. This is indeed the case. The results can be described 
rather well by the following quadratic expression: 

= 0.16 + 85712(0) ^ 0.16 + M^R^{0) . (88) 
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(Note that ^ MpR^ for A = 0.1). Five dots correspond to the hve values of the exponent 
/?: the largest action conies from the theory with the largest exponent, see Fig. 9. 

The total wormhole action (which is of main interest for us) also grows as Mpi?^(0), but with 
a greater coefficient: 

5totai~5 + 1.7M2i?2(0) . (89) 

Because of the computational difficulties, we have not performed the calculation for the values 
of exponent greater than 500. At /3 < 500 the logarithmic terms which are taken into account 
in S'totai (but not in S^) also give a considerable contribution. Therefore the quadratic £t for 
S'totai is much less accurate than the £t for S^, and the coefficient 1.7 in (|59D is obtained with 
a rather limited accuracy. However, the behavior of the solutions allows us to make a plausible 
assumption. We expect that with the further growth of the exponential factor the held /(r) will 
become frozen near /q, and the total wormhole action will approach our expression for the action 
(^), (1^) for the Giddings-Strominger-Lee wormhole, 

^ (l - ^ (l - ^) , (90) 

This expression is consistent with our approximate quadratic £t (P9|). 


5 Global symmetries and Planck scale physics 

5.1 Kaluza-Klein wormholes 

As we have seen, it is extremely difficult to increase the wormhole action by changing the effective 
potential of the scalar held. However, as a result of our investigation we have learned that in those 
cases when we were able to make the action large, the value of the action could be estimated by 
a simple expression S'totai ~ 1.34i?^(0)Mp. Remember also that the action becomes almost three 
times larger if one does not include the inner boundary contribution. This suggests that one can 
obtain a very large action if there exists some reason why the wormhole throat should be large. 
Indeed, as we have mentioned in the Introduction, we may not have any problems with axions if 
the ehective coupling constants of the operators violating the global symmetry are smaller than 
g-i90_ (|88D suggests that this happens if the radius of the wormhole throat becomes greater 

than 15Mp^. We were unable to make the throat that large by changing the ehective potential 
of the scalar held, but there exist other possibilities to do so. 


32 





Indeed, we have assumed that our space remains four-dimensional and that gravitational 
interactions are described by the standard Einstein theory at all length scales. Meanwhile, each 
of these assumptions may be wrong. 


First of all, according to Kaluza-Klein theories, the number of dimensions of space-time is 
much greater than 4, but space-time becomes effectively four-dimensional at i? > Rc-, where Rc is 
the radius of compactihcation. What if Rc S> Mp Then our equations should be considerably 
modihed at i? < Rc, which may lead to the wormhole throat of a large size i?(0) ~ Rc- 


It is not easy to test our hypothesis in any realistic theory, but we may play with a toy model. 


First of all, at i? S> our equations should coincide with our original equations (|^ , (^). 
Meanwhile we will assume that at R Rc space-time becomes ten-dimensional, like in string 
theory. In this case at R ^ Rc the charge conservation equation instead of n = 2tt‘^P9'R^ gives 
n = ^ i 2 R^^ ■ (W^e assumed that 6 dimensions form a sphere S'g of radius Rc and area and 


took into account that the area of a sphere Sq of radius R is Modihed equations equations 


), (1^ ) look as follows: 
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72n^R 


12 
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= 0 , 


= 0 . 


(91) 


(92) 


Note also that the factor 3 in front of in the hrst of these equations was replaced by 9, which 
corresponds to space-time with d = 10. 


The idea of a phenomenological description of possible wormhole solutions in this situation is 
to solve equations which at R ^ Rc look like (0), (|5B|), but at R Rc look like (PT|), (P^). It 


can be achieved, e.g., by introducing an interpolating factor (1 -|- which changes the 

9-dimensional volume, and, correspondingly, the conservation law, n = 2Tr‘^p9'R^(l + 

This equation gives n = 2 'k‘^P9'R? at i? S> Rc, and n = at i? -C Rc- After this and some 

other obvious modihcations the interpolating equations can be represented in the following way: 
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0, 

(93) 

0 . 

(94) 


We have solved these equations numerically for various values of Rc, see Fig. 10. As expected, 
the radius of the wormhole throat i?(0) was found to be approximately equal to the compactih- 
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cation radius Rc. Of course, our investigation of Kaluza-Klein wormholes cannot be considered 
conclusive. Still it indicates that the global symmetry breaking problem may disappear in the 
theories where the radius of compactihcation Rc is sufficiently large. 


5.2 One-loop effects in quantum gravity 


Another possibility is related to the TZ^ corrections which may appear in the effective Lagrangian 
or in equations of motion of the gravitational field. The simplest example is the conformal 
anomaly, which gives the contribution to the gravitational equations. 

Here 


— 2 71 + 27Z7Z^i, — -g^uTZ^ , 


(3) H = 

fliy 


7Z^ 7Z\y 


Inn 


liu 




(95) 


The parameters Hq and Mq are of the same order as Mp, but they can be much smaller than Mp 
if there are many matter fields (of spin 0, 1/2, and 1) contributing to the conformal anomaly. 


Equation for the scalar field (|^ remains unchanged with an account of conformal anomaly, 
but the gravitation field equation (^) acquires some new terms. To get an idea of the possible 
influence of quantum corrections on the structure of wormholes, let us assume for simplicity that 
Hq -C Mq, so that the second term in (|^ can be neglected. In this case the gravitational 
equation looks as follows: 


- 1 + 


SttR^ 

3M| 


V{f) + 
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Stt^PR^ 
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This equation on the throat yields 




{R'^ - 1 )" = 0 


(96) 


R^(n) = f/-2 - (v( n 4 - i 


(97) 


In the limit Hq S> Mp an additional term — does not alter our wormhole solutions. 

However, for smaller values of Hq the character of our solution changes dramatically, see Fig. 11. 
The throat of the wormhole becomes considerably wider, and the interval of r where i?' -C 1 
becomes very small. Finally, this interval disappears altogether, and for Hq < 2Mp regular 
wormhole solutions with = 0 cease to exist. In other words, even small quantum gravity 

corrections can lead to absence of wormhole solutions! 
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5.3 String-inspired models 


We have mentioned in Sect. that the hrst wormhole solution was in fact obtained in the version 
of the axion theory where instead of the pseudoscalar axion field one has the field In 

the same paper [^] Giddings and Strominger have obtained a family of wormhole solutions in the 
string-inspired version of the theory of the held with the effective action 

Sduai = I d^S{K - Ko) . (98) 

Here 0 is the dilaton held, 0 is a phenomenological parameter; /3 = 2 in string theory. For 
simplicity, we used here units in which = 1. (Wormhole solutions in a more general class 
of theories have been obtained later by Coule and Maeda p^). The solution for R{r) in this 
theory does not depend on (3. It coincides with the corresponding solution in the theory (^). 
However, the situation with the dilaton held is more complicated. The solution can be written 
in the following form: 

^ g-f0(o) cos arccos^^J^^ . (99) 

For Id < ^ this equation describes the wormhole solution with 0(r) which gradually increases 
at large r from its maximal value 0(0). However, for P > -^ this regime becomes impossible. 
Indeed, in this case 0(r) becomes inhnitely large at the point where arccos becomes 
equal to |. At this point derivatives of the scalar held diverge, and the action becomes inhnitely 
large. The conclusion of ref. pj was that in such a situation there are no regular wormhole 
solutions in this theory. 

Does it mean that the gravitational ehects cannot make the held massive in the theory 
(p8|)? We have already discussed a similar situation in the theory of a scalar held without 
symmetry breaking, where the total action was inhnite, and the conclusion was that the ehects 
responsible for the global symmetry violation are suppressed only by the part of the action coming 
from a small vicinity of the wormhole throat. Thus one could argue that the global symmetry 
violation could occur in the theory (^) as well, even despite the absence of the wormhole solutions, 
if one considers a small part of the conhguration (|9^ near the wormhole throat. 

However, the main reason why this argument could work for the theory of a scalar held was 
the existence of two vastly diherent length scales. The wormhole throat had a nearly Planckian 
size ~ whereas the typical scale on which the scalar held signihcantly changed was much 

greater, of the order of m~^. Therefore it was possible to pack many wormhole throats inside the 
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region of size m This is not the case for the conhgurations (^). In the most interesting case 
of the string theory with (3 = 2 the held ^(r) becomes inhnitely large at the point where 

^ 1.27/2(0) . (100) 

Thus the total size of our held conhguration almost coincides with the size of its throat. Since 
each of such conhgurations has inhnite action and is not of a wormhole type, we do not expect 
them to lead to global symmetry breaking. 

This does not necessarily imply that there are no wormhole solutions in string theory. The 
model considered above does not contain any potential V{<(>). Also, it is very hard to associate 
the value of the dilaton held 0, which is typically assumed to be of the order of Mp, with the 
parameter /o ~ 10^^ GeV. Nevertheless, this model clearly shows that the existence of wormhole 
solutions in the axion theory is by no means automatic. In this model the wormhole solutions 
disappear as soon as one considers ehects associated with the dilaton held. 


R{r) = 


R(0) 


COS 
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6 Topological suppression of wormhole effects in gravity 
and string theory 


All our previous results have been obtained by an investigation of particular solutions which may 
or may not appear in diherent theories. However, there is one general reason which may lead to 
suppression of wormhole ehects. These ehects lead to the change of topology of space by creating 
a universe capable of carrying a global charge away from our space. On the classical level such 
processes simply cannot occur. Our use of Euclidean methods to describe such processes still 
needs to be fully justihed. But even if these methods are valid, there is an easy way to suppress 
the probability of the wormhole formation in the Einstein theory. 

In order to explain it, we remind that the standard Lagrangian of the Yang-Mills held in 
QCD originally was written as , but later it was recognized that one can add to this 

Lagrangian the term ^i^F^yF^'^. This term does not modify the Yang-Mills equations, but it 
gives a contribution to the nonperturbative processes involving change of topology of the Yang- 
Mills held. 

A similar situation occurs in the Einstein theory of gravity, where the standard Einstein 
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Lagrangian --^IZ can be supplemented by two different topological terms, 

s = - j . (loi) 

Here *Tl^„xs = 


The last two terms do not give any contribution to equations of motion, and therefore do not 
change the theory at the classical level. Therefore, just as the term , these two terms 

can be considered as an integral part of general relativity. The first of these terms is very similar 
in its nature to the term . It contributes to the effective potential determining the 

value of the axion field. Fortunately, it is not expected to lead to any problems with strong CP 
violation. The effects induced by this term are related not to the wormhole physics, but to the 
Abelian Eguchi-Hanson instantons, and typically they are exponentially small, being suppressed 
hy e where e is the electromagnetic coupling constant . 


Meanwhile, the Gauss-Bonnet term gives a nonvanishing topological con¬ 

tribution to the wormhole action, 

^topoi = . (102) 

Therefore it may control the strength of the global symmetry breaking. 


This term was considered in the early works on wormholes @]. The constant 7 was called 
there a topological coupling constant. Just as the 6 parameter, the value of this constant in 
gravitational theory is not determined a priory. 


It is useful to remind the reason why the effects of the Yang-Mills instantons are suppressed 

1 h 

by the factor e dF. The semiclassical Yang-Mills action ■^FahF°‘° is not topological. Therefore 
the variation of this action produces equations of motion which have instanton solutions. The 
semiclassical action calculated on these solutions is known to suppress the instantons F^j^i, = 
\e^u\&Fx& as follows: 


^-s 


= e 



(103) 


where ly is the winding number of the gauge configuration. 


In gravity things seem to work differently, but the results are very similar. It is the Einstein 
action with matter which gives us the equations of motions. Those equations have wormhole solu¬ 
tions. For the Euclidean wormhole geometry (^) with the wormhole radius R(r) the topological 
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(Gauss-Bonnet) contribution to the action is given by the following integral: 


^topoi = + n^) = ^ I . (104) 


After the angular integration the integral becomes 

5'topoi = Y Jo ~ ■ (105) 

Keeping in mind that the Gauss-Bonnet part can be brought to a form where it is a total 
derivative, we can rewrite this integral as 

^topoi = Y Jo ^ “ ^i?'^(r)]^^oo = 7 • (106) 


In performing this calculation we did not use any particular form of the wormhole solution; we 
have used only the fact that -R'(O) = 0 and R'{oo) = 1. Thus, independently of any suppression 
calculated in the previous sections of this paper, there exists an additional exponential suppression 
of the wormhole-related effects in quantum gravity by the factor e"'*'. This agrees with the result 
obtained in p. 


It is important that this additional suppression is equally related to the total probability of 
the wormhole formation and to the values of the vertex operators. Indeed, for all wormhole 
solutions which are known to us the derivative R'{r) approaches its asymptotic regime R'{r) = 1 
at r comparable with the radius of the wormhole throat. Thus, the integral in eq. (|106|) rapidly 
converges to 7 in the vicinity of the wormhole throat. It means that all vertex operators become 
suppressed by the factor e~'^. 


Note that the value of parameter 7 is arbitrary; it does not change any experimentally tested 
predictions of general relativity. If one takes 7 > 190, all our problems with axions in quantum 
gravity immediately disappear. One may or may not like having a large parameter 7 in gravita¬ 
tional theory, but it is certainly not forbidden, and it solves the problem of the global symmetry 
breaking. 


Still it would be nice to hnd some reasons why this parameter should be large. One of the 
ideas is to consider string theory and to study an analogous topological term there. 

In string theory it is considered plausible that the Gauss-Bonnet term appears at the level of 
a' corrections in a specihc form since it is related to the Green-Schwarz mechanism of cancellation 
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of anomalies. One expects that this part of stringy corrections has the following form 

a! 


C-. 




..stringy - l 

Here a' is a fnnction of the fnndamental dilaton held if 


( 107 ) 


a = 


9 I 
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Onr normalization corresponds to the standard normalization of the Yang-Mills Lagrangian, 


.^YM — 


a 


IQk"^ 
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and we consider as the “rnnning” gange conpling constant - 57 -^. 

Qq 9 


What follows is a direct generalization of onr pnrely gravitational calcnlation. The only 
difference is the presence of the fnnction a'{x). Unfortnnately, we do not know mnch abont the 
dependence of a' on x, which makes the resnlts which we are going to obtain less rigorons bnt 
perhaps still rather plansible. First of all, since we stndy spherically symmetric conhgnrations, 
we snppose that a' depends only on r and does not depend on the angnlar variables. The integral 
becomes 


rcxD 

‘S'topol 127r / diT — r 

JO 9 [r) 


R'\r){l - R'\r)) . 


( 110 ) 


As we have seen in the previous sections, on a sufficiently large distance from the wormhole 
throat the wormhole geometry becomes undistinguishable from the geometry of a flat Euclidean 
space because its scale factor R{r) becomes almost exactly equal to r, and its derivative R'{r) 
rapidly approaches 1. Typically it happens at the distance of the same order of magnitude as 
i?(0); precise value of will not be particularly important for us. It is important, however, that 
at a sufficiently large r > the term (1 — R!‘^{r)) in the integral ( p.lO|) becomes very small, which 
implies that the total value of the integral is dominated by integration in a region r < r^. Note 
also, that on all our solutions we had (1 — R!‘^{r)) > 0 and i?" > 0. Therefore one can represent 
the integral (| 110 |) in the following way: 
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where is some average value of the gauge coupling constant in the region 0 < r < defined by 


eq. (| 111 |) . Since we expect to be of the same order of magnitude as i?( 0 ), and R{0) should be 
determined by typical stringy length scale (just as the natural scale for i?(0) in the Einstein theory 
was given by Mp^), we will identify with the typical value of the gauge coupling constant gstr 
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on the stringy scale. The subsequent evaluation of the integral in (|111| ) goes exactly as in (|106|) , 
because i?' 1 on the boundary r = r^. This gives 


5'topol — 


9stT 


( 112 ) 


Note that this is a topological contribution to the action, which practically does not depend on 
the detailed form of the wormhole solutions (this dependence is concentrated in our dehnition of 




5'str). 


This part of the action is a precise analog of the one-instanton action ^ in the Yang-Mills 
case. This leads to an additional suppression of the wormhole-induced effects by the factor 


p ‘S'topol 


= e 


(113) 


We would like to stress that whereas in the Yang-Mills theory the suppression of the instanton 
effects comes from the semiclassical action, the suppression of the wormhole effects described 
above did not came from the action of Einstein gravity with matter but from the topological 
term in the action which in string theory appears at the level of a' corrections. Thus, the term 
‘S'topoi = ^ appears here in addition to the usual wormhole action. 


Let us estimate the numerical value of Stopoi in the realistic theories. We do not really know 
the value of the effective coupling constant on the wormhole throat (i.e. approximately on 
the stringy scale). The simplest idea would be to identify with the gauge coupling constant 


related to the grand unihcation in supersymmetric GUTs, ~ ^ |^. This would 

give S'topoi 163. Thus, the topological suppression alone can be strong enough to eliminate 
the effects of the rank Eve operators (Remember that the coupling constant g^ for 


/o ~ 10^^ GeV should be smaller than 10“®^ ~ HI]-) To get the factor required 

to suppress the most dangerous term giMp{^ -|- <h*) we need a slightly smaller gauge coupling 
constant, 

1 

< — . (114) 

30 ^ ^ 


Air 


Note that the factor ^ is completely analogous to the topological coupling constant 7 which 

3str 

we discussed in the case of pure gravity. In that case 7 could take any possible value, and it 
was not quite clear whether it is natural or not to take it as large as 190. In the case of the 
string theory the condition 7 > 190 corresponds to a very natural constraint ^ This 

g2 

requirement seems quite reasonable since the effective gauge coupling constant ~ T can 
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slightly decrease on the way from the GUT scale 2 x 10^® GeV to the stringy scale, which is much 

2 

higher. Moreover, one can solve all problems even with ^ ~ ^ if in addition to the topological 
action m) one takes into account the standard contribution to the action S which we have 
studied in the previous sections. 


It is very hard to discuss these issues in the absence of a well established string-inspired 
phenomenological theory. Nevertheless we will try to make some simple estimates. With this 
purpose we should remember relation between the stringy gauge coupling constant Qstr, the stringy 
mass scale Mstr (i-e. the mass of the hrst massive stringy excitation), the stringy length scale Zgtr 
(compactihcation radius), and the parameter a' in the heterotic string theory |l3l|, |3^ : 

2 _ 2 _ gstr Mp 

^ o! Zstr 


Mstr — 


(115) 


Here a' is the effective value of the parameter a' on the scale of the wormhole throat. Therefore 
the topological action ^ can be also written as 


‘S'topol - 


8M 

9stT 


^a'Ml = = 71 



(116) 


Let us hrst concentrate on the expression S'topoi = tt • Uncertainty in the value of 

translates into the uncertainty of M^tr- If one simply takes Mstr ~ obtains Mstr ~ 

1.66 X 10^® GeV, which then leads to the estimate for the stringy unihcation scale Ugtr ~ 4 x 10^^ 
GeV. However, this implies the existence of a large gap between the stringy unihcation scale and 
the SUSY GUT unihcation scale 2 x 10^® GeV. Therefore there is a tendency to assume that for 
some reason the stringy scale is in fact considerably smaller than Mstr 1.66 X 10^^ GeV m- 


In order to obtain a sufficiently strong suppression due to topological ehects it would be enough 
to have Mgtr ^ 1-5 x 10^® GeV, which is quite consistent with the present ideas about the value 
of Mstr- 


Moreover, this constraint on Mstr can be somewhat relaxed and reduced to Mstr ^ 2 x 10^® 
GeV if one takes into account the standard (nontopological) contribution to the wormhole action. 
Indeed, as we have seen, the usual contribution to the GSL wormhole action can be represented 
as 1.34Mpi?(0)^ (1^. This result was quite consistent with our results for the theory with the 
exponential potential and with simple dimensional estimates suggesting that the action of the 
wormhole with the radius of the wormhole throat i?(0) should be of the order M|i?(0)^. Note 
that the topological contribution also has the same structure (although with a slightly smaller 
coefficient), S'topoi = f^p^str- If fhc wormhole solutions exist at all in the string theory, one 
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may expect that the wormhole throat -R(O) should be greater than the “elementary length scale” 

4tr- This suggests that the nontopological part of the action ~ Mpi?(0)^ should be greater than 

If this is correct, the total action of the wormholes, including the topological part, should 

be about two times greater than the topological contribution. In fact, if one simply adds to 

the topological action the GSL action 1.34Mpi?(0)^ for R{0) ~ Igtr, the value of the 

topological action almost triples. In such a situation one can expect that the total action should 

become greater than 400 even if one takes ^ ~ “ 4 ^ ~ This is more than snfficient to 

solve the problem of the global symmetry violation. One would get the total action S > 190 for 
2 , 

^ which looks like a very safe bet. 

Thus, instead of the absolutely incredible hne-tnning of the valnes of possible coupling con¬ 
stants of the interaction terms breaking the global symmetry in the axion theory, onr estimates 
gave us rather mild constraints on the gange coupling constant and on the stringy mass scale 
Mstr- This provides a natnral possibility to make stringy gravity compatible with the existence 
of the light axion in Nature. 
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7 Discussion 


One of the main obstacles on the way of development of quantum gravity is the problem of 
experimental verihcation of its predictions. Indeed, gravitational interactions between elementary 
particles become strong only at energies comparable to the Planck mass, Mp ~ 1.2 X 10^9 GeV. 
Such energies are far beyond our reach. 

Fortunately, there exist some indirect ways to test quantum gravity experimentally. For 
example, it is quite possible that a consistent theory of gravity requires supersymmetry. Then one 
may study different versions of the theory by investigation of the properties of light superpartners 
of the graviton. Another possibility is related to cosmology, which provides us with experimental 
data originated at the very early stages of the evolution of the Universe. 

In addition to that, there may exist some nonperturbative gravitational effects which may 
have important experimentally testable consequences even at very low energies. Some of these 
effects have been investigated in this paper. We have found, under the assumptions specihed 
in this paper, that gravitational effects strongly violate global symmetries in a wide class of 
theories, including the theory of axions. Surprisingly enough, this strong violation occurs even if 
one drastically modihes the effective potential of the theory, for example, if one multiplies it by 
exp where the factor C can be as big as 10^. 

Of course, it is quite possible that our methods based upon Euclidean approach to quantum 
gravity are inadequate. However, we must admit that when we began our investigation we 
expected that it will be very easy to fix this problem either by using a formulation where the 
global symmetries become local, or by finding a simple modification of the theory which leads to 
wormholes with a very large action. We have found that it is almost impossible to do so in the 
context of the standard Einstein theory in a four-dimensional space. 

However, as soon as we allowed ourselves to modify the theory of gravity or the properties 
of space at the scale about lOMpT^ we have found many different possibilities to improve the 
situation. One of them is the possibility that our space is compactified, with a compactification 
radius Tc > We have seen also that wormholes may simply disappear if one takes into 

account conformal anomaly, or if one considers certain string-inspired axion models. In addition 
to all these effects, we have found that in string theory there exists a specihc strong suppression 

of topology change by the factor e ®str = e . This is a topological effect which does not 

depend on many particular details of wormhole conhgurations. Our estimates show that with an 
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account taken of topological terms, the effects related to the global symmetry breaking become 
strongly suppressed if the string mass scale is sufficiently small, Mstr ^ 2 x 10^® GeV. Equivalently, 
the problem of the wormhole-induced global symmetry breaking disappears if the gauge coupling 
constant is sufficiently small on the stringy scale of energies, ^ < jg- These values of parameters 
are quite consistent with the existing picture of stringy phenomenology. Therefore at present we 
do not see any reason to reject the theories with global symmetries. 


On the other hand, we have found that the existence of (approximate) global symmetries 
and the possibility to solve the strong CP violation problem by the Peccei-Quinn mechanism are 
very sensitive to the choice of the theory of quantum gravity. If the axions in the mass range of 
nia ~ 10“® eV will be discovered experimentally |^, it may give us important information about 
the structure of space and properties of particle interactions at the Planckian scale. 
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Appendix: Numerical Methods of Finding Wormhole So¬ 
lutions 


Finding the wormhole solutions numerically was a rather complicated problem. We have found 
that the behavior of the solutions was extremely sensitive to the choice of initial conditions. 
Whereas it was relatively easy to hnd solutions with small action S, in the most interesting cases 
where the action was large the standard numerical methods failed. Therefore it was necessary 
to develop a more advanced method of calculations. There is a chance that our method can be 
useful for solving other problems as well. Therefore in this Appendix we will hrst describe the 


standard numerical methods of obtaining solutions to the differential equations (|73|) , (|74[) , and 
then we will describe our method. 


Standard Shooting Method 


The hrst method we try is the standard shooting method |^. We hx some value for F’(O). Then 
we hnd A(0) using ([55|). This involves solving a cubic equation for A^(0). 


A^(0) - A®(0)1/(F(0)) 
If there are no roots, then eqs. 


F2(0) 


= 0 


(117) 


, (|74|) do not have a solution for the hxed value of A(0). Eq. 
( p.l7|) has at most two roots. We hnd that the bigger root does not lead to a solution satisfying 
boundary conditions at inhnity. Hence we pick the smallest root. 


Once we have the initial conditions we use any numerical ordinary diherential equation solver 
(for example Runge-Kutta method). According to the behavior of the solution at inhnity we 
modify our guess of F{0). Using a simple iterative procedure we can determine the correct value 
for F(0). 

First we have used Runge-Kutta method to obtain the numerical solution once the initial 
conditions (boundary conditions at x = 0) are set. However, for solutions with large action, 
the solution has to be computed to a very high accuracy. We have found that the Runge-Kutta 
method does not work for actions above 6. Bulirsch-Stoer method, which we used next, works 
for actions below 10. As a result, simple shooting method is inadequate for our purposes. The 
reason it does not work is that K(0) has to be determined to a very high accuracy, otherwise 
the numerical solution will not come close to the boundary conditions at inhnity (at x ~ |R). 
On a computer, K(0) can be determined only to machine precision. The difference between 
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F(0) rounded to machine precision and the correct -F(O) causes an exponentially large deviation 
between the numerical solution and the correct solution for large x, hence making the standard 
shooting method fail. 

Improved Shooting Method 

With the standard shooting method we cannot get a numerical solution to closely approximate 
the correct solution for the wormholes with large action. The idea of an alternative method is to 
sacrihce the correctness of the numerical solution at some x, and use the gained freedom to make 
it closer to the correct solution. The disadvantage is that the resulting numerical solution does 
not satisfy the differential equation at all points. The advantage is that it is close to the correct 
solution. 

We solve the differential equation in several stages. The hrst stage uses the standard shooting 
method with the Bulirsch-Stoer method as the ordinary differential equation solver. After A(0) 
has been determined within machine precision, the hrst stage is completed. Every time a solution 
is generated with the Bulirsch-Stoer method, the program remembers at which x the solution 
rapidly veered up or down so that it is clear it will not satisfy the boundary conditions at inhnity. 
Once the stage is completed, the program recalls solutions F^p{x) and Tdown(a^) which are the 
curves which have passed next to the correct solution for the longest time. Fy^p{x) has passed 
above the correct solution (so it veered up) and Tdown(a^) passed below the correct solution. 

The program Ends the maximum x for which the two solutions are still very close to each 
other: |Fup(a;) — Tdown(a^)| < e. Let this point be Xm- Then we have found a good approximation 
to the correct solution up to x = Xm- The second stage is to try to find an approximation to 
the solution for x > Xm- For the second stage, and all further stages, the initial conditions are 
different from the ones at the hrst stage. As in the case of the hrst stage, we are free to choose 
F{xm) with the restriction Fkp(a;m) > F{xm) > Fdown(a^m)- Notice that by choosing F{xm) we 
introduce a point in the numerical solution which does not satisfy the diherential equation. The 
error introduced is very small, though, since F-ap{xm) is very close to Fdown(a^m)- 

We also need to know A{xm)- We hnd it by linear interpolation between Aup(a;m) and 
Adown(a^m) once F{xm) is set. In the same way we determine F\xm)- As already mentioned, 
the jump in F{x), A[x), F\x) at x = Xm is small and can be made as small as we like by reducing 
the parameter e. During the second stage F{xm) is determined to machine precision by the same 
iterative technique used in the hrst stage to determine E(0). Then we proceed to the third stage. 
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and so on, until we reach a; ~ after which the solution is nearly constant: F[x) ~ Fq, x oo. 

The method described above can be interpreted in a way similar to the interpretation of the 
standard shooting method. We aim and shoot for the boundary conditions at inhnity. Having 
missed several times, we follow the trajectory the best arrows followed and approach the target a 
little closer. Then we shoot from the new position. After hnally hitting the target we present the 
set of trajectories of different arrows as one nearly correct trajectory which a single arrow could 
have followed all the way from the initial position to the target. Being a bad shot, we would not 
have been able to shoot the arrow to follow the correct trajectory all the way from the boundary 
at a; = 0 to the boundary at inhnity. However, using the method described above we can hnd a 
good approximation to the correct arrow trajectory. 
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Figure Captions 


Figure 1: (a) The geometry of the Giddings-Strominger-Lee wormhole with initial surface 
Sj which is and hnal surface Sj which is x S^. This wormhole may describe the tunneling 
from R^ to R^ x S^. (b) Extended solution which connects two asymptotically Euclidean regions 
r —±cxD. 

Figure 2: The distribution of logj^gF as a function of log^gp for the wormhole solution in the 
theory with an effective potential — /q)^. Here p = rMpy^, F = 

Figure 3: Wormhole geometry for the potential ^{p—foY- Here p = rMp A = RM-p^j^. 

Figure 4: Total action as a function of — InFo = In ^ + | In ^ for A = 0.1. The black dots 
and the line show the total action with the boundary term at the throat, the grey ones give 
the corresponding values for the action without the boundary term. Note that S'totai ~ 15 for 
/o ~ 10^2 GeV. 

Figure 5: Plot of log^g ^ 1°^ different values of the exponent and /o = 10^® GeV . 

Figure 6: Scale factor A{p) for different values of the exponent and /o = 10^® GeV. 

Figure 7: Derivative of the scale factor A{p) for different values of the exponent and /o = 10^® 
GeV. Wormhole geometry approaches geometry of a flat Euclidean space when A'{p) approaches 
1 . 


Figure 8: The total action in the theory with the exponential potential V^f) = {p — pY 

as a function of the exponent j3. Black dots give the action with the boundary term, grey one 
give the action without the boundary term at the wormhole throat. 

Figure 9: The action obtained by integration in the region near the wormhole throat (where 
H'(p) = R!{r) < 0.9) for different values of the exponent j3. 

Figure 10: Behavior of the radius R of the wormhole solution at small r in our model of a 
10-dimensional Kaluza-Klein theory with different values of the compactification radius Rc. All 
quantities are given in the Planck units Mp = 1. 

Figure 11: Behavior of the radius A{p) of the wormhole solution at small p in the model 
taking into account conformal anomaly. Note that regular wormhole solutions with A'(0) = 0 
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(i.e. with -RhO) = 0) disappear for Hq < 2Mp, even though in this regime the coefficient in 
front of the conformal anomaly term is still very small. 


52 



